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Abstract 


A point matrix kernel for radiation transport, developed by the 
transmission matrix method, has been used to develop buildup factors 
and energy spectra through slab layers of different materials for a 
point isotropic source. Combinations of lead-water slabs were chosen 
for examples because of the extreme differences in shielding properties 
of these two materials. 


Introduction 


The transmission matrix method as applied to 
radiation transport tHeory was developed by Yarmush 
et.al. (ref. 1). This development was applicable 
^or"uncharged particles from an infinite plane 
source in slab geometry. Angular distributions were 
approximated by angular cosine polynomial expansions 
and energy distributions were approximated by a 
group structure. 


Description of the Point Kernel 


As a result of the development (ref. 1) , the 
transmission through a plane slab is 

4 pl - t(x>s a) 

where : 

T(x) = the transmission matrix 

S = a source vector representing the 

source angular and energy distribu- 
tion 

^ =s a response vector representing the 
transmitted angular and energy dis- 
tribution. 

Recent work (ref. 2) has shown that one can obtain 
a point kernel from equation (1) . This kernel is 

K “ — ^ [T(x) - xT'(x)]. (2) 

p c a™ 2 

After application of the transmission matrix 
for shields consisting of two different materials, 
one has (ref. 1) 

$ = T 2 (x 2 )[I - R 1 (x 1 )R 2 (x 2 )]' 1 T 1 (x 1 )S 

(3) 


where : 

R i (x i ) 


* reflection matrix for material i 


a * the fraction of total thickness in 
1 material i. 


By definition material one is located next to the 
source. In a manner similar to analyses which lead 
to equation (2), one can write, for the transmission 
matrix through two slabs, 

4ttx 2 4> pt 31 t 1 ~ x I)*|t 2 (x 2 )[I - R 1 (x 1 )R 2 (x 2 >] 

XT i (x i ) }-‘ ' (4) 

Performing the differentiation results in 

4ttx 2 4 = [ t 2 ( x 2 ) M(x 1’ x 2^1 (x 1 ) 

- x 2 T'(x 2 )M(x r x 2 )T 1 (x 1 ) 

/ 

- x 1 T 2 (x 2 )M(x 1 , x 2 )T 1 (x 1 ) 

- T 2 (x 2 )M( Xl , x 2 )[x 1 R{(x 1 )R 2 (x 2 ) 

+ x 2 R 1 (x 1 )R 2 (x 2 )]M(x 1 , x 2 )T 1 (x 1 )] S . ( 5 ) 

where 

M.( Xl ,x 2 ) - [I - R 1 (x 1 )R 2 (x 2 )] l . 

At first indication this expression appears awkward 
for practical application; however, since the 
response effects from a double ref lection for gamma 
rays can usually be neglected, one can write 

4TTx 2 <£ pt = [T 2 (x 2 )T 1 (x 1 ) - x 2 T'(* 2 )T 1 (x 1 ) 

- x 1 T 2 (x 2 )T 1 / (x 1 )] S. (6) 

Equation ( 5 ) can readily be modified to include any 
number of slabs of materials. Neglecting the 
reflection effects, one has 

4TTxV = (I - x ^D[T n (x n ) ••• T 2 (x 2 )T 1 (x 1 )] S. 

P (7) 

Note that these expressions are matrix expressions, 
and the commutation property with respect to matrix 
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multiplication does not hold in general. 

The transmission matrix is given as (ref. 1) 

T(x) = 4C + _1 e' Ax D(x) (8) 

where 

D(x) * (B + + B e -Ax C C + " 1 e"^)" 1 

= (I - B + - 1 B_e- Ax C -C+ - 1 e- Ax + . • 

The matrices B + , B_, C + , and C_ are functions of 
the material scattering and absorption properties 
only and not of material thickness. These opera- 
tors are directly related to the radiation trans- 
port equation. Therefore the thickness occurs only 
in the exponential matrix which is a diagonal 
matrix. In most practical problems the double ex- 
ponential and higher products can be neglected and 
equation (8) becomes the asymptotic form 

T .(x> = 4C/ 1 e' Ax B + _1 (9) 

and the first derivative of T(x) becomes 

T^(x) = - 4C + 1 Ae _Ax b” 1 . (10) 

Therefore equation (2) becomes 

K pt = ~^~2 C + 1(1 + Ax)e _Ax B -1 . (11) 

TTx " r 

By similar analyses equation (6) becomes 


K-edge of the photoelectric cross section for lead. 
The effect is observable only for very thin sheets 
of lead. 

Figure 7 illustrates the energy fluence build- 
up factors for a combination of three slabs of lead 
and water. These results show the effect of moving 
a lead slab through a water shield. Figure 8 shows 
the energy spectra for a two mean free path lead 
slab. This figure illustrates how effectively the 
lead slab removes the low energy photon build-up 
in the preceding water slab. 


Conclusion 

The point matrix kernel has been shown to be a 
relatively simple and effective means for calcula- 
ting energy spectra and build-up factors for com- 
binations of materials in slab geometry for a point 
isotropic source. This particular kernel should 
become an effective tool for parametric and optimi- 
zation studies . 

In this work only gamma ray shielding is con- 
sidered; however, there is no fundamental reason why 
the procedure will not also work for any particles 
whose behavior can be described by a linear trans- 
port equation. The only distinction would be in 
the calculation of the particle cross sections. 
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Similar results can be obtained for n slabs. 


Results for Multi-Layered Slabs 


Results for layered combinations of lead and 
water are illustrated in figures 1 to 8. The choice 
of materials was not made for any particular shield 
configuration but rather to illustrate a combination 
of shielding material extremes. Therefore they 
should represent a good test for a shielding calcu- 
lation. In all cases a source of 1 MeV is used. 

Figure 1 shows the energy fluence build-up 
factor for combinations of lead slabs followed by 
water slabs. Figures 2 and 3 are the differential 
energy spectra for total shield thicknesses of five 
and ten mean free path lengths. 

Figure 4 shows the energy fluence build-up for 
combinations of water slabs followed by lead slabs. 
As the water slabs become thicker one can note the 
sharp drop in the build-up factor due to even very 
thin slabs of lead. The shape of these curves are 
similar to those of Kalos (ref. 3). This effect is 
further illustrated in the energy spectra in 
figures 5 and 6. These results show the effects of 
the photoelectric cross section of lead. In figure 
6 one can note the "filtering" effect due to the 
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ENERGY FUUENCE BUILDUP FACTOR 




Energy fluence build-up factors for lead-water 3. 
slabs . 


Energy spectra for lead-water slabs of ten mean 
free paths total thickness. 




Energy spectra for lead-water slabs of five 4. Energy fluence build-up factors for water-lead 

mean free paths total thickness. 








Energy spectra for water-lead slabs of five 
mean free paths total thickness. 



Energy spectra for water-lead slabs of ten mean 
free paths total thickness. 


7. Energy fluence build-up factors for water-lead- 
water slabs. 



Energy spectra for water-lead-water slabs of 
ten mean free paths total thickness for a con- 
stant thickness of two mean free paths of lead. 






